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We study bridging transitions between spherically and cylindrically shaped particles (colloids) of 
radius R separated by a distance H that are dissolved in a bulk fluid (solvent). Using macroscopics, 
microscopic density functional theory and finite-size scaling theory we study the location and order 
of the bridging transition and also the stability of the liquid bridges which determines spinodal lines. 

The location of the bridging transitions is similar for cylinders and spheres, so that for example, at 
bulk coexistence the distance Hb at which a transition between bridged and unbridged configurations 
occurs, is proportional to the colloid radius R. However all other aspects, and, in particular, the 
stability of liquid bridges, are very different in the two systems. Thus, for cylinders the bridging 
transition is typically strongly first-order, while for spheres it may be first-order, critical or rounded 
as determined by a critical radius R c . The influence of thick wetting films and fluctuation effects 
beyond mean-field are also discussed in depth. 


I. INTRODUCTION 

It is well known that effective pair interactions between dissolved particles may be induced by the presence of the 
solvent. A much studied example of this phenomenon is a solvent-mediated depletion interaction between colloids 
which occurs when their size is substantially larger than that of the solvent particles pQ. In this case, the entropy 
gain due to excluded volume give rise to an effective short-ranged attraction between the colloids determined by the 
diameter of the solvent. The force of depletion may then lead to clustering of the solute particles and plays a crucial 
role in the stability of colloidal and colloidal-polymer systems. 

Another type of effective force may occur when the solvent forms a liquid or gas bridge between the colloidal 
particles. The details of this depend sensitively on the wetting properties of the colloidal particles, the distance 
between them and also how close the solvent is to its own fluid-fluid phase boundary. While such bridging induced 
interactions require finer thermodynamic tuning than the force of depletion, they are stronger and also longer-ranged 
since they are determined by the size of the colloidal as opposed to solvent particle. Bridging between colloidal 
spheres was originally studied using macroscopic arguments mu and simple effective interfacial models mm- In such 
macro/mesoscopic models the transition between unbridged and bridged configurations must be first-order owing to 
the different symmetries of the interfacial profiles. However, more microscopic approaches based on density functional 
theory (DFT) have shown that the scenario is more complex [7HID]. In particular, as shown recently 10j. the 
phenomenon is much more similar to capillary-condensation and criticality in parallel-plate geometries. This means 
that for spherical colloids the formation of a liquid bridge only occurs via a first-order transition for sufficiently large 
sphere radius R > R c . Surprisingly for R < R c the formation of a bridge is continuous and occurs via the smooth 
merging of the wetting layers surrounding each colloid - something which is beyond the reach of mesoscopic models. 

In this paper we compare and contrast the character of the bridging transition between spheres and between parallel 
aligned cylinders to which much less attention has been devoted m- We show that there are major quantitative 
differences between these two types of bridging which are apparent at macroscopic level, within microscopic DFT and 
beyond mean-field when fluctuation effects are considered. Of particular importance is the stability of liquid bridges, 
which are very different for the spheres and cylinders, both at and off bulk coexistence. Thus for example at bulk 
coexistence for cylinders our DFT studies indicate that, at least away from the immediate vicinity of the bulk critical 
point, and in contrast to the case of spheres, there is no critical radius R c , so that a first-order bridging transition is 
always present. We show that at two phase coexistence, the shapes of the minimal surfaces forming the bridge lead 
to universal results for the location of the bridging transition for large particle sizes. This is closely analogous to the 
Kelvin equation for condensation in capillary slit geometries and implies, in the macroscopic limit, that the distance 
H between the particles at which bridging occurs scales linearly with the particle size R. Corrections to this resulting 
from the presence of thick complete wetting layers are also predicted and tested using DFT, as is the case when the 
solvent fluid is off two-phase coexistence for which a non-linear phase boundary and spinodal are constructed. The 
finite-size rounding of the bridging transition due to fluctuation effects is also discussed. 
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FIG. 1: Macroscopic phase diagram for first-order bridging transitions between completely wet cylinders immersed in a solvent 
at two-phase coexistence. The line represents the bridging coexistence curve Hb — ttR. The bridged phase remains metastable 
for all H > Hb. 


II. MACROSCOPIC THEORY 

In this section, we formulate a simple macroscopic theory for bridging transitions between two identical parallel 
cylinders and between a pair of identical spheres. The distance between these, measured centre-to-centre (for spheres) 
and axis-to-axis (for cylinders), is H while the radius in each case is R. These are immersed in a vapour at temperature 
T < T c with T c the critical temperature and chemical potential /i which may be at or away from bulk two phase 
coexistence /i = /i sat (T). For the most part we will assume that the spheres/cylinders are completely wet by liquid 
so that the contact angle for each is zero although it is straightforward to generalise the analysis to partial wetting 
which is also discussed briefly. If the spheres/cylinders are far apart they are each coated with a layer of liquid 
the thickness of which is determined by R and the range of the intermolecular forces. However, the forces and the 
thickness of the wetting layers need not be specified in the macroscopic regime. For smaller separations H on the 
other hand it is possible that liquid condenses in the region between the spheres/cylinders forming a bridge between 
them. This reduces the area exposed by the wall-gas interface but increases that of the wall-liquid and liquid-gas 
interface. A free-energy balance between these two configurations occurs at H = Hb(R), determined by the radius 
R, corresponding to the macroscopic bridging transition. At this level the transition is always first-order since the 
configurations are always distinct. However, the line of the bridging coexistence curves Hb(R) is different for the 
cylinders and the spheres because the shapes of the liquid bridges, which are minimal surfaces at bulk coexistence, 
have very different properties. In the following, let q w i, q wg and qi g represent surface tensions between wall-liquid, 
wall-gas and liquid-gas interfaces, respectively, which are assumed to be the same as those of the planar surfaces. 
The differences between the free energies of the unbridged and bridged configurations arise simply from the surface 
tensions and the appropriate area exposed. We consider this for cylinders and spheres separately. 


A. Bridging between cylinders 

Let us begin by considering the case of cylindrical colloids that are completely wet (0 = 0) and a solvent that is at 
bulk coexistence. In the unbridged state the excess free-energy per unit length of the cylinder is F = 47rq wg R coming 
from two separate cylinders. In the bridged configuration the cylinders are connected by a liquid film consisting of two 
parallel flat liquid-gas interfaces. The curvature of these is trivially zero since we are at bulk coexistence. These flat 
interfaces meet the cylinders tangentially since the contact angle is zero. The excess free-energy for this configuration 
is F = 2ttR^ w i + 2irR^ wg + 2H^\ g , implying that the bridging transition occurs when 


Hb = 7 tR . 
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FIG. 2: (Color online) Macroscopic phase diagram for first-order bridging transitions between completely wet cylinders 

immersed in a solvent off coexistence. The solid line represents the bridging coexistence curve H^/L — 2R/Lf(R/L) whilst the 

dashed line is the spinodal curve H sp i n /L = 2w ^ at which the bridging him pinches. 



This is equivalent to the Kelvin equation for condensation in capillary slits but now written in terms of the cylinder 
separation rather than the undersaturation. Note that this simple linear relationship is to be expected from dimen¬ 
sional considerations since R and H are the only macroscopic lengths in the problem. The macroscopic phase diagram 
for bridging transitions occurring at bulk coexistence is therefore particularly simple and is shown in Fig. 1. For 
fixed R and H > the liquid bridge is metastable with respect to the unbridged configuration, however there is no 
spinodal determining the limit of stability since flat interfaces may always span the two cylinders. 

The result is readily generalised to the case of partial wetting for which the contact angle 6 > 0. The liquid-gas 
interfaces of the bridge phase are still fiat but now meet the cylinders at the contact angle rather than tangentially. 
A simple calculation yields 


—^ = ( 7 r — 20) cos 0 + 2 sin 0 , ( 2 ) 

R 

which is symmetric about 0 = 7 r/ 2 . Notice that H 5 = 2 R for neutral wetting, 0 = 7 r/ 2 , which means that the 
macroscopic bridging transition no longer exists as the cylinders are in contact. The value of H 5 is maximum for the 
case of complete wetting (0 = 0) or complete drying (0 = n). For the latter scenario, a bridge of gas forms between 
the cylinders when it is immersed in a solvent of bulk liquid. 

For the case of complete wetting one can also consider the influence of the singular part of the wall-gas surface 
tension arising from complete wetting layers. For example, for the case of short-range intermolecular forces, the 
wall-gas surface tension is modified to m 


7wg — 7wi A 7ig A 


In 

R M 



(3) 


where £* is the bulk correlation length of the solvent in the liquid phase. Making use of this result, the bridging 
transition is shifted to 


IT 6 «7ri2 + 7r&ln(|) , (4) 

implying that the wetting films enhance the onset of the transition by increasing the effective cylinders size. Similar 
corrections occur for dispersion-like forces where they are algebraic of order Rz [IMS]. The same corrections occur 
for complete drying films. Thus for systems with short-ranged forces, the expression for H 5 is of the form @ but 
with the bulk gas correlation length £ g replacing 

Finally let us return to the simple macroscopic theory for complete wetting and consider what happens off bulk 
coexistence when the chemical potential or pressure departure from two-phase coexistence Sp = p sa t — p is non-zero. 
In this case the liquid-gas interfaces forming the liquid bridging film are not flat but rather arcs of circles of radius 
L = 7 i g /Sp as required by Laplace’s equation. The free-energy balance argument must now also include a term from 
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FIG. 3: Illustration of the length-scales characterizing a macroscopic bridged configuration for spheres, see text for details. 


the cross-sectional area of the bridge which is proportional to Sp and a simple calculation shows that the transition 
occurs when 


H b (R,L)=2Rf(j/j , (5) 

where f{x) is a scaling function satisfying the implicit equation 

sin_1 = (TW (f (2 + ^ “ A/( 1 + *) 2 -* 2 / 2 ) . (6) 

It is easy to show that /(0) = 7t/2 while /(oo) = 1 so that in the limit R/L —>• 0 and R/L —>> oc one obtains the 
linear results H b « 7 ri? and H 5 « 2i?, respectively. The first limit therefore recovers our earlier result for bridging 
transitions at bulk coexistence. The second limit means that when the colloid radius is larger than the Laplace radius, 
the bridging transition occurs very close to the contact distance. Interestingly, this generalised phase diagram now 
includes a spinodal curve, see dashed line in Fig. 2, H s = 2 \/R? + 2 LR, which marks the limits of stability of the 
bridging phase, corresponding to the p inching of the bridge film in which the two arc circles of radius L meet. Notice 
that the spinodal curve H 5 « y/SLR for small R/L so that the stability of the bridge increases as coexistence is 
approached. However in the opposite limit, R/L —oc the curves H s and H b have the same asymptote so that the 
stability of the bridge is reduced. This means that when R^> L, the bridging transition occurs near contact and is 
weakly first-order. This is the macroscopic mechanism for the disappearance of the bridging transition as T T c , 
since in that limit, L —>• 0 due to vanishing of the surface tension 7 i g . 


B. Bridging between spheres 


Next consider the bridging transition between two spherical colloids, each of a radius R , immersed in a bulk vapour 
at /i = /i sat . Due to the last condition, the mean curvature of the macroscopic shape of the liquid bridge must be 
zero, i.e., must form a minimal surface. It is an easy task from a variational principle to show that he minimal surface 
of a revolution is a catenoid formed by rotating the catenary y = a cosh ^ about the revolution x-axis connecting the 
centres of the spheres. The bridging transition occurs when there is a balance of surface free energies of the bridged 
and unbridged states balance equivalent to the condition 

Scat _ 7wg ~ 7wl /y\ 

Swi 7ig 

where S cat is the area of the catenoid and S w \ is the total area of the spheres covered by liquid in the bridged 
state. Following the illustration in Fig. [3j the catenoid connects with the spheres at x = ±h and therefore S w \ = 
4:7rR(R + h — Hb/2). The catenoid has area 


Sc. 


f h 2 
= 2tt a / cosh 

Jo 


fx\ 

a ( 2 h\ 

( — ) dx = 2ira 

- sinh ( — + h 

\aJ 

.2 V a ) 


( 8 ) 


and substitution gives 


a 

R(2R + 2h-H b ) 




= cos 0 


(9) 
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where we have used Young’s equation. The parameters h and a are determined from the contact conditions 


a cosh 


s/R 2 - {h - F b /2) 2 


( 10 ) 


and 

. _i [ H b /2 - h 
tan — ; = 

y?? 2 - (/i - i? 6 /2) 2 

where in the last equation we have allowed that the liquid-gas interface meets the spheres at an arbitrary contact 
angle. These equations have the solution 


— tan 1 

sinh (-) 


V«/J 


■e, 


(ii) 


f = «(«), (12) 

where the function a(6) has to be determined numerically and is shown in Fig. [4] The separation H 5 is largest for the 
case of complete wetting (or complete drying), for which <a(0) « 2.32 in an agreement with the analysis of [7b The 
separation H 5 is smallest for neutral wetting (0 = 7 r/ 2 ) for which again = 2 R meaning that no bridging occurs. 
The linear relationship between H 5 and R for Sp = 0 is similar to that found for cylinders. However, the stability of the 
liquid bridge is completely different. In contrast to cylinders, the metastable region occurring for H > is of finite 
extent and the spinodal line H sp [ n (R) lies close to the coexistence curve H^{R). The spinodal occurs not because the 
bridges pinch, as in the case of cylinders, but because for large H the catenoid has no solution touching the spheres. 
This is illustrated in Fig. [ 5 ] where we show the coexistence curve H^R) ~ 2.32 R and spinodal H sp [ n (R) « 2.38 R for 
the case of complete wetting. In the inset we show the dependence of the length-scale a characterising the catenoid 
as a function of the sphere separation H. Notice that the value of a/R = 0.363 at the spinodal is close to that at the 
bridging transition for which a/R = 0.476. More generally the spinodal line is always a linear function of the radius, 
^spin(^) = Oi sp [ n (6)R, but has a coefficient that depends on the contact angle as shown in Fig. [4j This shows that 
the extent of the metastable region is largest for the case of 6 = 0 and smallest for 0 = n/2. None of this behaviour 
is present for cylinders, where the spinodal line is only important off bulk coexistence. The close proximity of the 
spinodal line for spheres is suggestive that the macroscopic bridging phase diagram at Sp = 0 may be more strongly 
influenced by microscopic effects, than for the phase diagram for cylinders. 

Finally we mention that similar to the cylindrical case, there are non-analytic corrections when complete wetting 
layers coat the spheres. For systems with short-ranged forces these lead to 


Hb 

R 


c*(0) + 




(13) 


in which the pre-factor of 2 in the log correction arises from the curvature of the sphere. This means that the 
corrections due to complete wetting films are smaller for spheres than for cylinders. Once again, for complete drying 
the bulk correlation function £/ is replaced by £ g . 


III. DENSITY FUNCTIONAL THEORY 

Within classical density functional theory m the equilibrium density profile p(r) is found by minimising the grand 
potential functional 

= F[p] + j dr Vi(r) - lA p( r) • ( 14 ) 

where Vi( r), i = s,c, is the external field due to the pair of spherical or cylindrical colloids, respectively, and p is the 
chemical potential. It is convenient to split the intrinsic free energy functional T into the contribution due to the 
ideal gas and the remaining excess part 

F[p\ = ‘Fid \p\ + F'exM , (15) 

where [p] = ksT J drp(r) [ln(A 3 p(r)) — l] with A the thermal de Broglie wavelength that can be set to unity. 

Following the standard perturbative approach, the excess part of the free energy can be separated into repulsive 
and attractive contributions: 


•^exM = ^hs[p] + \ J drp(r) J dr'p(r')«a(|r - r'|), 


(16) 
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FIG. 4: (Color online) Contact angle dependence of the parameter a(0) appearing in the macroscopic prediction Hb = olR for 
spheres (blue) and cylinders (red, uppermost) in which case the parameter a is given by Eq. No bridging is possible for 
6 — The dashed line represents the contact angle dependence of <a S pm(60 determining the spinodal H s p i n = ot sp in(R) for 
spheres which lies close to the coexistence line. 



FIG. 5: (Color online) Macroscopic phase diagram for first-order bridging transitions between completely wet spheres immersed 
in a solvent at bulk coexistence. The solid line represents the bridging coexistence curve Hb ~ 2.32 R whilst the dashed line is 
the spinodal curve H sp i n /L « 2.38 R. The inset shows the dependence of the parameter a, characterizing the catenoid shape 
a cosh as a function of H. The termination of the curve determines the spinodal at which the catenoid no longer touches 
the spheres. 


where u a (r) is the attractive part of the intermolecular potential. In our calculations we take this to be that of a 
square-well (SW) potential 

usw = ^hs + u a , (17) 

where uus is the hard sphere potential with diameter a and 

( 0; r <c r, 

Wfl(r)a< -6] a<r<Xa , (18) 

I 0; r > \a , 


where we set A = 1.5, as is usual. The hard-sphere potential contribution to the free energy functional is included 
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FIG. 6: Schematic behavior of the solvation force showing possible scenarios for bridging transition between spheres or cylinders: 
(a) first-order bridging, (b) critical bridging and (c) rounded bridging. 


within Rosenfeld’s fundamental measure theory m 

^hs [p\ = ksT J dr<F({n a }). 


(19) 


The free energy density <£> is a function of a set of three independent weighted densities {n a } that can be expressed 
as two-dimensional integrals mm- Here we adopt the original Rosenfeld prescription, although other forms of <£> 
are available m- 

In our tests of the macroscopic theory we will focus only on the case of complete drying (0 = n) for spherical and 
cylindrical colloids in a bath of saturated liquid. This is most easily modelled by using a purely hard-wall repulsion 
for the spheres and cylinders which induces drying at all temperatures T <T C . 

The minimisation of Eq. (14) leads to the Euler-Lagrange equation 


SF 

k B T In A 3 p(x, r) + — + V s (x, r) = p , 

dp 

for spherical colloids (with r the perpendicular distance to the x axis) and 

SF 

k B T In A 3 p(x, z) + —— + V c (x, z) = p , 
dp 


( 20 ) 


( 21 ) 


for cylindrical colloids. 

The equations (J20| and (21) are solved iteratively on a two-dimensional grid with a uniform spacing of 0.05cr in all 
directions. We adopt a Cartesian coordinate system for cylinders (exploiting a translation invariancy in the direction 
along the cylinders) and a cylindrical coordinate system for spheres (taking advantage of the axial symmetry). In 
both cases, an overall finite-size box is of a rectangular shape of dimensions L x x L z with L x chosen up to 500 a and 
L z up to 200 a. By determining the equilibrium density profile, the excess (over ideal) grand potential U ex can be 
calculated. The effective potential between the pair of immersed particles a distance H apart is given by 


W(H) = n ex (H) - 2QW , 


( 22 ) 


where flix is the excess grand potential corresponding to a system where only one colloidal particle is present. The 
derivative f s = — dW/dH defines an effective force of solvation between the spheres or cylinders. Using this microscopic 
formalism we can go beyond macroscopics and classify the order of bridging transitions. As the sphere/cylinder 
separation is reduced a first-order bridging may occur at H = i^, where the value H^ = H^R) is a function of 
the sphere/cylinder radius. For large R we anticipate that this H^R) is a near linear function of R in accord with 
macroscopic expectations. At the separation H^ the potential W(H) has a kink meaning that the solvation force 
f s (H ) is discontinuous and exhibits a finite jump at H = H 5 . It is, however, conceivable that the line of first-order 
bridging transitions H^R) terminates at a critical point R = R c at which the solvation force is infinitely steep, i.e. 
df s /dH = 00 . The bridging transition is then critical. Finally, for R < R c the solvation force f s (H ) is a continuous 
function from which one can identify a rounded bridging transition at H = where the curvature of the solvation 
force changes sign as given by d 2 f s /dH 2 = 0. The possible nature of bridging transitions are illustrated in Fig. [ 6 j 
Our hard-wall colloids, which are completely dry, are immersed in saturated liquid (solvent) with the temperature 
set to T = 0.89 T c where /c^T c = 0.9 e is the bulk critical temperature. For each value of the colloidal separation 
H and different radii R we minimize the grand potential functional U[p(r)] on a two-dimensional grid to obtain the 



















FIG. 7: (Color online) Phase diagram for cylindrical bridging transitions obtained for T — 0.89 T c at a two-phase coexistence. 
The symbols represent DFT results obtained from a numerical solution of Eq. (21). The red dashed line is a straight line with a 
slope of 7r in accord with the macroscopic prediction given by Eq. 0 - The black dotted line is a fit containing a log correction 
due to complete wetting by gas around the cylinders parameterized results as Hb = ttR + a\nR/a. We find with a — 3.34 a in 
good agreement with the predicted value a = 7r^ g « 3.52. 



FIG. 8: (Color online) Two-dimensional density profiles of the saturated liquid near to two infinite cylinders each of a radius 
R c = 70 a whose centers are separated by a distance H = 230 a. The density profiles correspond to a first-order bridging 
transition at which the unbridged state (a) and the bridged state (b) coexist. 


equilibrium density profile. This is done by starting from two different initial configurations corresponding to a bridged 
and an unbridged state where the bridge is formed by a single film of a low-density gas connecting the colloids. In 
standard fashion, convergence to distinct configurations then indicates the presence of stable and metastable states 
and the crossing of stable branches in the Q-H diagram determines for first-order bridging transitions. The 
absence of distinct branches indicate that the formation of a bridge as H is continuous corresponding to either critical 
or rounded transitions. 

In Fig. [7] we show the bridging phase diagram for the values of vs R obtained for cylinders with radii ranging 
from R = 3a to R = 70a. There is excellent agreement with the macroscopic prediction H^ = ttR for R > 10 a and 
even better agreement over the entire range of R values when we include a logarithmic correction as expected from 
the influence of the complete drying layers. It is surprising that the latter expression provides almost a perfect fit to 
our DFT data down to radii of a few of a units, despite the mesoscopic nature of the prediction. Fig. [8] shows two 
coexisting density profiles corresponding to bridged and unbridged states for R = 70 a. Notice that the liquid-gas 
interfaces of the gas bridge connecting the cylinders are near planar in agreement with macroscopic expectations. 
A remarkable feature is that for the present temperature a first-order bridging transition occurs down even to the 
smallest radius considered. We have also determined bridging phase boundary off bulk coexistence, results of which 
are shown in Fig. [9j In contrast to the case of bulk coexistence, the H^-R dependence is no more linear but rather 
follows the macroscopic law given by Eq. (J5| with two linear regimes appropriate for small and large R. The bridging 
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transition remains first-order down to about R c « 3cr, but exhibits continuous (rounded) transition below. In Fig 
we show two density profiles corresponding to the phase boundary for R = 5 cr (first-order bridging, bridged state} 
and R = a (rounded bridging). 

For spheres, the phase diagram is similar but the role of a critical radius R c is more pronounced presumably because 
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FIG. 9: Phase diagram for cylindrical bridging transitions obtained for T = 0.89 T c off coexistence showing comparison between 
DFT (symbols) and macroscopic theory (solid line) as given by Eqs. § and <§• Here, the Laplace radius L « 5 a. The bridging 
transition remains hrst-order down to microscopic radii but exhibits a critical point below the radius R = 3a (see the text). 
Note that in contrast to the case of coexistence shown in Fig. [7] where the macroscopic theory was confirmed by analyzing the 
slope of the phase boundary for sufficiently large R , here the full comparison is made for the whole range of radii considered. 
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FIG. 10: (Color online) Two-dimensional density profiles of the supersaturated liquid near to two cylindrical colloids each of a 
radius R = a (a) and R = 5 a (b) corresponding to the bridging transition. For R = 5 a the transition is first-order and the 
boundary of the bridging film is near circular with a radius of curvature L « 5 a in line with the macroscopic considerations. 
For R — cj the transition from a bridged to an unbridged state is continuous. 



of the weakened stability of the bridging film. This is illustrated in Fig. [Tl] showing the phase diagram obtained for 
the same temperature T = 0.89 T c with radii ranging now from R = 3 a to R = 25cr. In this case we find that a 
first-order bridging transition only occurs for R > R c with R c ~ 10 a. In the first-order regime there is once again 
near perfect agreement with the macroscopic prediction H^ = 2.32 R. Logarithmic deviations from this are not at all 
significant as the radius R is large. Also recall that the anticipated pre-factor of any log correction is only half of 
that predicted for cylinders. Coexisting density profiles and the fit to a catenoid are shown in Fig. [12]. For R < i? c , 
as the separation H is reduced, a gaseous bridge forms continuously between the spheres without the occurrence of 
any singular behavior in the effective potential. Thus the solvation force f s (H) shows no jump and is a continuous 
function. However, we may identify a rounded bridging transition for R < R c which occurs when there is a point of 
inflection f' s '(H ) = 0, as shown in Fig. [I3| together with the solvation force behaviour for R > R c exhibiting a finite 
jump. Remarkably, the macroscopic prediction H 5 = 2.32 R extends through the critical point at R c and predicts 
rather accurately the location of the rounded transition (cf. Fig. 11). 
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FIG. 11: Phase diagram for spherical bridging transitions obtained for T = 0.89 T c . The symbols refer to the DFT results: the 
filled squares denote points for which the bridging transition is first-order, while the empty squares correspond to a continuous 
bridging transition. The line represents a linear fit with slope 2.38 close to the predicted macroscopic value 2.32. The critical 
radius is estimated as R c ~ 10 a. 



FIG. 12: (Color online) Two-dimensional density profiles of the saturated liquid near to two spherical colloids each of a radius 
R = 25 a whose centers are separated by a distance H = 64 a. The density profiles correspond to a first-order bridging transition 
at which the unbridged state (a) and the bridged state (b) coexist. A fit to a catenoid shape is shown (black line). 


IV. DISCUSSION 

In summary, we have used macroscopic theory and microscopic DFT to study the bridging transition between 
cylinders and between spheres. The transitions are both qualitatively similar to condensation occurring in capillary- 
slits and pores but are quantitatively different to each other due to the very different stability properties of the 
liquid bridge. In particular for spheres, DFT studies show that there is always a critical radius R c below which the 
bridging is rounded while for cylinders, the limit of stability depends on the Laplace radius of the bridging interface, 
so that the transition remains first-order at bulk coexistence and away from the critical temperature. This echoes 
the stability properties of the liquid bridge within the macroscopic theory where it is only for spheres that, at two 
phase coexistence, the spinodal and coexistence lines lie close to each other. For bridging between spheres our DFT 
studies indicate that the critical radius R c increases with temperature. For example repeating our numerical analysis 
we have found that at T = 0.94 T c the value of the critical radius for spheres is R c ~ 30 cr. Indeed, for a near critical 
temperature T = 0.98 T c we could not detect first-order bridging transition at all for any values of the colloids radii up 
to R = 70 a. This is consistent with expectations of finite-size scaling that R c should increase with £/ or £ g . Indeed, 
for spheres recall that even at the macroscopic level the distance between the bridging transition and its spinodal at 
bulk coexistence, is extremely small with, H sp i n — « 0.0 6R. Minor corrections to either of the lines, H^R) and/or 

iLspin(^), coming from microscopic length-scales, may cause these to cross, which would indicate the location of a 
critical radius. A natural candidate for such a microscopic length is the wetting layer thickness which is oc £/ In R/&, 
as predicted by Eq. pi). Assuming that H sp [ n (R c ) — Hb(R c ) oc £/ In R c /£i yields a value for R c which scales with the 
bulk correlation length but with a large pre-factor. 

The role of a critical radius R c is at least quantitatively different for the case of cylinders. In this case we are 
able to find first-order bridging transitions that persist much closer to T c and down to small radii. This is consistent 
with the above finite-size scaling argument since at bulk coexistence there is no macroscopic spinodal line. However, 
discussions concerning the existence of a critical radius R c and its behaviour as T both for spheres and cylinders 
increases are somewhat academic, since beyond the mean-field, fluctuation effects destroy all true phase equilibria, so 
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FIG. 13: (Color online) DFT results for the solvation force f s (H) of a fluid with short range forces adsorbed between two spheres 
as a function of separation for two different radii, (a) R = 25a > R c , first-order bridging (blue line), (b) R — 7.5 a < R c , 
rounded bridging (red line, rightmost). The dotted vertical line represents the macroscopic prediction for Hb/R — 2.32. 



FIG. 14: (Color online) Schematic phase diagram for bridging transitions and typical configurations between completely wet 
(dry) cylinders taking into account the effect of fluctuations. The dotted line represents the location of the macroscopic phase 
boundary Hb = ttR whilst the dashed line represents the cross-over region of width ~ e - 2 ^^ig R centered around Hb(r) over 
which the bridging transition is rounded. In the cross-over region we have highlighted the typical size of the domains of bridged 
and unbridged regions. 


that all bridging transitions are rounded. For bridging between spheres the volume available to the gas or liquid bridge 
is finite, meaning that the first-order boundary Hb(R) must be smeared over a region proportional to the volume of 
the system [21,, i.e. of order R~ 3 . Similarly for cylinders, where the bridge extends along the cylinders, we anticipate 
that the rounding is sharper and is over order e ~ 27r ^s R which is Boltzmann’s factor corresponding to the probability 
to form a liquid-gas domain, since the geometry is pseudo one-dimensional. In this case, near to the macroscopic line 
Hb = ttR the bridge breaks up into domains of bridged and unbridged states of extension £ « e 2?r ^7ig^ ? similar to 
the pseudo-symmetry breaking in 2D Ising strips, see Fig. [l4j It is this rounding which provides the true mechanism, 
beyond mean-field, for the disappearance of even rounded bridging transitions as the temperature is increased to T c . 

In this work we restricted our microscopic (DFT) study for the models in which the fluid particles interact purely 
repulsively with the walls (cylinders or spheres). Although from a macroscopic perspective the situation in which the 
fluid interacts attractively with the walls and a bulk phase is vapour rather than gas is completely analogous, one 
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expects a richer phase behaviour of such systems with prewetting and layering phenomena also involved when a more 
microscopic treatment is used. Moreover, the high-density phase near the walls will exhibit much stronger excluded 
volume effects (not included within the macro- or mesoscopic approaches) so that the comparison between analytic 
predictions and numerical DFT results is expected to reveal somewhat larger deviations. 

It would also be interesting to consider other geometries such as crossed or non-parallel cylinders for which the 
shape of the bridging film is more complicated requiring a three dimensional DFT analysis. In this case we anticipate 
that the bridge may undergo further phase transitions and associated interfacial fluctuation effects. An example of 
this would be reducing the opening angle between two slightly misaligned, near parallel, cylinders, in which case the 
length of the bridging film must diverge. This is similar, but not identical to filling transition for fluids adsorbed in 
wedge geometries. 
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